Survival probability of a diffusing particle constrained by two moving, absorbing 

boundaries 
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We calculate the exact asymptotic survival probability, Q, of a one-dimensional Brownian particle, 
initially located at the point x G (— L, L), in the presence of two moving absorbing boundaries 
located at ±(L + ci). The result is Q(y, A) = ^2^ = _ 00 {—i) n cosh(ny) exp(— n 2 A), where y = cx/D, 
X = cL/D and D is the diffusion constant of the particle. The results may be extended to the 
case where the absorbing boundaries have different speeds. As an application, we compute the 
asymptotic survival probability for the trapping reaction A + B — » B, for evanescent traps with a 
long decay time. 
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Physical systems described by partial differential equa- 
tions with moving boundary conditions are ubiquitous in 
nature [l| . Unfortunately, such equations are notoriously 
difficult to solve. The case of a single moving boundary 
is often amenable to analysis, for example by tranform- 
ing to the moving frame, but the case of more than one 
moving boundary is, in general, intractable. 

First-passage problems are another field of research for 
which there are relatively few exact results @, Q • The 
simplest such problem, for which some exact results are 
available, is that of a Brownian particle (i.e. random 
walker) moving in the presence of one or more absorbing 
boundaries Q. The case of a single boundary, moving 
at constant speed, can be solved exactly 0] but, to our 
knowledge, the survival probability of a single Brownian 
walker in the presence of two moving boundaries with 
different velocities had not been solved up to now. 

In the present paper we apply "backward Fokker 
Planck" methods to solve this problem exactly. As an 
application, we consider the one-dimensional trapping re- 
action A + B — > B, where the density of traps, p, decays 
exponentially with time, and obtain the exact asymptotic 
form of the final A-particle density in the limit where the 
decay-time, r, of the traps is large. 
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FIG. 1: Brownian walker, starting from the origin, in a 
linearly-expanding cage. 

We consider a Brownian walker moving according to 
the Langevin equation X(t) = rj(t), with initial condition 
X(0) = x, where ij(t) is Gaussian white noise with mean 
zero and correlator (r)(ti)r)(t2)) = 2D6(t\ — £2)- The 



particle is flanked by two receding, absorbing walls with 
coordinates ±(L + ct), as shown in Fig. [I] We can derive 
a backward Fokker-Planck equation for the probability, 
Q(x,L,t), that the particle has survived up to time t 
having started at position x G (—L,L). We note that, 
after infinitesimal time At, the particle will have have 
moved to position x + AX, and the walls will have moved 
to positions ±(L + cAt). It follows that Q(x,L,t) = 
(Q(x+AX, L+cAt, t—At)), where the average is over the 
distribution of the spatial increment AX. Expanding to 
first order in At, using (AX) = and ((AX) 2 ) = 2D At, 
yields the backward Fokker-Planck equation 
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The infinite-time result is obtained by setting the time 
derivative to zero. It is convenient to introduce the vari- 
ables y = cx/D and A = cL/D to represent the dimen- 
sionless initial positions of the particle and the walls. 
Eq. ([T]) then reads, at infinite time, 



d 2 Q 

dy 2 



^ = 
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where —A < y < A, subject to the absorbing boundary 
conditions, Q(±A, A) = 0, and Q(y, 00) = 1, since if the 
particle starts at one of walls it is immediately absorbed, 
while if the walls are initially infinitely far away the parti- 
cle will survive with probability 1. A solution of Eq. @ 
satisfying the boundary conditions may be deduced by 
inspection, noting the symmetry of the problem under 
reflection, y — > — y; 



Q(y,X)= J2 (-l) n cosh(™,)e- 



(3) 
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Despite the simplicity of its derivation, Eq. ([3]) is, 
to our knowledge, a new result. It should be noted 
that this result seems difficult to obtain using conven- 
tional ("forward") Fokker-Planck methods. The back- 
ward Fokker-Planck method has eliminated the time- 
dependent boundary conditions from the problem at the 
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cost of introducing the initial wall parameter, L, as an 
additional independent variable. 

In Fig. [2] we present the results of numerically evaluat- 
ing the sum in Eq. ([3]) for various values of the starting 
coordinate, y. 
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FIG. 2: Survival probability Q{\) for various values of the 
initial coordinate y. 

For a particle starting at the origin, the survival prob- 
ability Q(0, A), written simply as Q(A), is 



Q(A) = J2 (- 1 )" c_ 



(4) 



To leading order for large A we may approximate the sum 
by the n = and n = ±1 terms: 



Q(A) ~ 1 - 2e" 



(5) 



This result agrees, to first order in e _A , with the ap- 
proximate infinite-time result, Q(X) ~ exp(— 2e~ x ), ob- 
tained using the method of Krapivsky and Redner [4J for 
a rapidly expanding cage. In this approach one writes 
the joint probabability distribution that the particle sur- 
vives till time i, and is located at x, in the approx- 
imate form 0| P(x,t) = [Q{t)/y/4irDt] exp(-x 2 /4Dt), 
i.e. one multiplies the free diffusion propagator by the 
probability, Q(t), for the particle to survive till time t, 
ignoring the boundary conditions at the walls, where the 
density is anyway small. The rate of change of Q(t) 
is minus the total probability flux through the walls, 
dQ/dt = 2D(dP/dx) x= L+ct- Integrating the resulting 
equation from t ~ to t = oo, using the method of 
steepest decents (valid for A ^> 1), gives the quoted re- 
sult. We see that this method only gives the leading 
departure from unity correctly. 

For small A, the leading-order behaviour can be ob- 
tained by rewriting Eq. (fj[]) using the Poisson sum for- 
mula > E5T=-ao f( n ) = Efcl-oc K^k), where f(k) is the 



Fourier transform of the the function f(n). This gives 

I — oo 

0(A) = a/t E e-^ 2 *" 1 ) 2 ^ (6) 



fc=-c 



The leading behavior at small A is 



Q(A) 



X 7r/A exp(-7r 2 /4A), 



(7) 



which contains an essential singularity at A = 0. We are 
unaware of any approximate methods for which even the 
leading behavior, Eq. ((7]), can be recovered. 

We have performed the summ in Eq. ^ numerically 
and plot the result together with the large and small A 
forms, Eq. (0) and Eq. ([7]), in Fig. H We see that the 
asymptotic forms describe the data well over a consid- 
erable range of A. This is readily understood on noting 
that the leading corrections to the asymptotic forms are 
of order exp(— 4A) and exp(— 97r 2 /4A) for large and small 
A respectively. 
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FIG. 3: Survival probability for a particle starting at the mid- 
point of a linearly-expanding cage with the asymptotic forms 
for small and large values of A, the dimensionless measure of 
the initial locations of the walls. 

To discuss the general case of walls with different 
speeds, we first consider the problem of a Brownian 
walker with drift ac following a path X(t) 1 X(Q) = x, 
with an absorbing boundary at the origin, and a receding 
absorbing boundary at position L + ct, where x € (0, L). 
The path X(t) of the walker satifies the Langevin equa- 
tion X(t) = r)(t) + ac, where rj(t) is again Gaussian white 
noise. We can write down a Backward Fokker-Planck 
equation for the survival probability of the particle at in- 
finite time which, making the same changes of variable 
as those used in Eq. (|2|), takes the form 



d 2 Q | JQ | dQ _ Q 
dy 2 dy d\ 



(8) 



subject to the boundary conditions Q(y = 0, A) = Q(y = 
A, A) — 0. In the limit A — > oo we impose as a boundary 
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condition the solution of Eq. (JH) with dQ/dX — 0, i.e. 
Q(y,X — > oo) = 1 — exp(— ay), which is the survival 
probability of a diffusing particle that drifts with mean 
velocity ac away from a fixed absorbing boundary or, 
equivalently, the survival probability of a particle with 
no drift in the presence of one wall receding at speed ac 
(see, for example, ref. Q and references therein). 

A solution to Eq. (JSj) satisfying all the boundary con- 
ditions is 



Q{y,\) = 



«y 



c -(n+a)y \ g -n(n+a)A 



(9) 



On shifting to the left both the path of the particle and 
the coordinates of the barriers by a distance (L/2 + act), 
we see that Eq. (j9|) solves the expanding cage problem for 
a particle with no drift, and with asymmetrically reced- 
ing walls — the left wall having coordinates — (L/2 + act), 
and the right (L/2 + (1 — a)ct), and may be shown to be 
equivalent to Eq. ([3]) when a = 1/2 (noting that the sepa- 
ration of the walls and drift velocities need to be doubled 
to show an exact mapping between the two problems). 

As an application of these results, we consider the trap- 
ping reaction, A + B — > B [;5|], with evanescent traps. We 
find the infinite-time survival probability of a particle, 
A, surrounded by an infinite sea of Poisson-distributed, 
evanescent traps, B, with initial density p, undergoing 
the trapping reaction in one dimension. By Poisson- 
distributed we mean that, at t = 0, each infinitesimal 
interval dx contains a trap with probability pdx. In par- 
ticular, the probability to find no traps in an interval of 
length L is cxp(— pL). The traps are evanescent in the 
sense that they randomly and independently disappear 
from the system in such a way that the overall trap den- 
sity decreases in a prescribed fashion, p(t) = f(t)p(0). 
The particle and traps both perform Brownian motion, 
with diffusion coefficients Da and D B respectively. This 
problem has recently been studied in detail for various 
functions f(t) for the case of a fixed target (Da = 0) 
Q. It was shown that, in this case, there is a non-zero 
infinite-time survival probability whenever f(t) falls off 
more rapidly that i -1 / 2 for large-t. Here we address 
the general and more difficult problem of a moving A- 
particle, and specialize to the case of exponentially de- 
caying trap density, f(t) = exp(— t/r), anticipating a 
non-zero infinite-time survival probability. This survival 
probability for a single A-particle also gives the fraction 
of A particles that survive if the initial state contains a 
macroscopic number of them. 

We approach the problem in the spirit of Bray and 
Blythe (7| , who considered the trapping reaction without 
trap decay, by finding upper and lower bounds on the 
survival probability and showing that they asymptoti- 
cally agree. In the calculation of both bounds we use the 
formalism introduced by Bray, Majumdar and Blythe 0. 
They define a quantity p(t), the mean number of differ- 
ent traps that would meet the ^4-particle up to time t , 



for a given ^-particle trajectory z(t), in a fictitious model 
where the A and B particles do not react. It satisfies the 
integral equation 

p(t)= [ dt' fi(t')G(z(t),t\z(t'),t'), (10) 
Jo 

where p(t) is the trap density, and G (z(t),t\z(t'),t r ) = 
exp [-(z(t) - z(t')) 2 /AD B (t - t')] /y/4nD B (t-t') is 
the trap diffusion propagator. Note that p[z] is implic- 
itly a functional of the A-particle trajectory, and that 
p(t = 0) = 0. The survival probability of the A particle 
is then given by [1, @| 



Q(t) 



(11) 



with the average taken over the A-particle trajectories 
with the usual Wiener measure. When the fraction of 
surviving traps is f(t), Eq. (fTO]) is modified thus [||: 



P 



dt' G(z(t),t\z(t'),t') 



o 



/(*') 



(12) 



where we now use p to refer to the initial trap density. 
This minor modification means that we may simply de- 
fine a new quantity <j>(t) = (i(i)/ f(t) in Eq. (fT2|) . solve as 
in the non-decaying case, and find that 



p(t) 



dt' 0(f')/(O- 



(13) 



We first derive an exact upper bound on the survival 
probability, as in 0, 0] , by arguing that the particle will 
survive longest if it remains stationary at the origin (the 
'target problem' (lpj). This is the so-called 'Pascal prin- 
ciple' [uj- We solve Eq. ([12]) with z(t) = 0, 



dt' 



<P(t') 



y/4arD B (t-try 



(14) 



and get the solution found in that is <ft — 

^Dst/ir) 1 / 2 which, using Eq. (|13|) with t = oo and 
f(t') = cxp(— t' /t), and Eq. (TTT]) . gives an upper bound 
on the eventual survival probability as 



Q < exp[-2(p 2 D B T 



,1/21 



(15) 



We can prove that this is a rigorous upper bound follow- 
ing the procedure outlined in [8|. If we write tfi = 4>q + 4>\ 
in Eq. ([T2"]). where (Ao is the solution to Eq. (TH| , we can 
show that 0i > 0, proving that <fi > <po. 

In [?[ , Bray and Blythe bound the survival probability 
from below by considering a notional box centred on the 
origin, from where the A-particle's trajectory is chosen to 
begin. They then select a subset of surviving trajectories 
by imposing three independent conditions that together 
guarantee the A-particle's survival: (i) The A particle 
does not leave the box up to time t; (ii) no trap enters 
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the box up to time t; and (iii) no traps were inside the box 
at time t = 0. These conditions undercount the number 
of possible surviving trajectories of the ^4-particle - for 
example, a trap may enter the box without trapping the 
particle - so the probability of fulfilling them underesti- 
mates the actual survival probability. We follow a similar 
line of argument here, but with a modification to adjust 
for the exponentially decaying traps: We allow the walls 
of the box to recede linearly. Then we may use the result 
obtained in Eq. (J4j) to determine the probability of sat- 
isfying condition (i) above. We consider, therefore, the 
survival at infinite time. The probability corresponding 
to condition (iii) follows simply from the Poisson prop- 
erty of the trap density, and for a box of length 2L is 
simply exp(— 2pL). 

To obtain the probability that no traps have entered 
the box, we consider each side independently and square 
the result. This is then equivalent to the survival proba- 
bility of a ballistic A particle in an infinite sea of traps, 
which was solved in [121 ]. The problem may be solved 
using a modified form of Eq. (fT2|) to allow for traps on 
only one side of the wall, a modification introduced in Q . 
With the linear trajectory z(t) — ct one obtains 



1 + erf 



•t exp 
= / dt'j>(t') 



y/4nD B (t-t') 

(16) 

We solve this for <j>(t) by Laplace transform and, using 
Eq. (fl3|) . obtain the probability 



exp -2m/ D b t 



(Da/ Db) 3 , which are satisfied for large enough r for any 
values of p, Da, and Db- 

Comparing the two bounds, (fl~5|) and (fl"9jl . we see that 
they agree, for large r, to leading order in the exponent. 
They therefore determine the exact large r asymptotics 
in the form Q ~ exp[—2(p 2 D bt) 1 ^ 2 + •■•], where the el- 
lipsis indicates subdominant terms. Note that the lead- 
ing term is independent of Da, just as was found for the 
time-dependent asymptotics of the survival probability 
for non-decaying traps Q. 

The fact that the bounds pinch indicates that the 
choice of linearly receding walls to obtain the lower bound 
is optimal for the case of exponentially decaying traps 
with a long decay time. Other forms for the trap de- 
cay function f(t) will require a different choice for the 
wall motion to optimize the bound. We recall that for 
non-decaying traps the bound is optimal for static walls 
0- 

In summary, we have obtained exact results for the 
infinite-time survival probability of a diffusing particle 
flanked by two receding, absorbing boundaries, including 
the case where the boundaries move at different speeds. 
We have used the results to compute the survival proba- 
bility of a particle diffusing in a sea of diffusing, evanes- 
cent traps in the limit where the trap decay time is large. 
Extensions to systems with general spatial dimensional- 
ity are possible and will be discussed elsewhere. 

The work of RS was supported by EPSRC (UK). 
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(17) 



that no traps have ever entered the box. 

We now combine these three factors to obtain a lower 
bound on Q: 



Q>exp -2p\ D b t 



per — 2pL — 



* 2 Da 
AcL 



where the last term comes from using the small-A result 
([7]) and we have neglected a logarithmic correction in the 
exponent coming from the pre-exponential factor in Eq. 
([HI). The result is valid for cL/Da <C 1 (corresponding 
to small A in the rescaled coordinates of Eqs. ([2])-([7])). 
This bound contains two free parameters, L and c, so we 
obtain the best lower bound by maximizing Eq. (|18|) with 
respect to both of these quantities, giving 



> exp -2( P 2 D b t) 



' :-pn 2 p 2 D A T)^ 



(19) 



This is valid as long as p 2 Dat ^> 1. 
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